History-Induced Critical Behavior in Disordered Systems 



John H. Carpenter, Karin A. Dahnien, Andrea C. Mills, and Michael B. Weissman 
University of Illinois at Urbana- Champaign, Department of Physics, 1110 West Green Street, Urbana, IL 61801. 



o 
o 

(N 



C/3 



Andreas Berger and Olav Hellwig 
Hitachi Global Storage Technologies, San Jose Research Genter, 650 Harry Road, E3 San Jose, GA 95120. 

(Dated: February 2, 2008) 

Barkhausen noise as found in magnets is studied both with and without the presence of long- 
range (LR) demagnetizing fields using the non-equilibrium, zero-temperature random-field Ising 
model. Two distinct subloop behaviors arise and are shown to be in qualitative agreement with 
experiments on thin film magnets and soft ferromagnets. With LR fields present subloops resemble 
a self-organized critical system, while their absence results in subloops that reflect the critical point 
seen in the saturation loop as the system disorder is changed. In the former case, power law 
distributions of noise are found in subloops, while in the latter case history-induced critical scaling 
is studied in avalanche size distributions, spin-flip correlation functions, and finite-size scaling of the 
second moments of the size distributions. Results are presented for simulations of over 10* spins. 

PACS numbers: 75.60.Ej, 75.70.Ak, 64.60.Ht, 75.60.Ch 
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Hysteresis occurs when a system far from equilibrium 
is driven by an external force. The system state then 
depends on the history of the system. In many sys- 
tems such as ferromagnets superconductors 13, 
and martensites 0], the response to the driving force is 
not continuous but occurs in discrete jumps which is of- 
ten referred to as "crackling noise", or specifically for 
magnets, Barkhausen noise. In magnets, broad ranges 
of power-law scaling of noise have been attributed to ei- 
ther a disorder induced critical point, such as found in 
the non-equilibrium zero temperature random-field Ising 
model (RFIM) 5|, or to self-organized criticality (SOC) 
13, as found in soft magnets. 

The disorder induced critical point found in the satura- 
tion magnetization curve, traced out as an external mag- 
netic field drives the system from one saturated state to 
the oppositely magnetized state and back, has been stud- 
ied in detail for the RFIM dHIl B and its existence 
has been experimentally confirmed [lOj. However, tak- 
ing magnets to saturation often is impractical due to the 
large magnetic fields required, so the behavior of subloops 
is of great interest to experiments and applications alike. 
The RFIM may be used to model subloops, and magneti- 
zation curves have been computed exactly in one dimen- 
sion JJJ and on a Bethe lattice fl3| and have also been 
collapsed near the demagnetized state using Rayleigh's 
Law Also, the idea of history acting as a source of 
effective disorder was recently introduced 14y . In this let- 
ter we report how the presence or lack of long range (LR) 
demagnetizing fields produces two distinct behaviors in 
subloops of the RFIM. In particular, in the presence of 
sufficiently large LR fields, subloops are found to resem- 
ble a SOC system. However, in the absence of LR fields, 
we find a critical subloop inside the saturation loop where 
the system history acts as a tuning parameter instead of 
the system disorder. Correspondingly, a modified scaling 



picture for subloops, consistent with the saturation loop, 
is introduced. For both cases we report new experimental 
results for subloops showing qualitative agreement with 
the two behaviors of the RFIM. 

In the non-equilibrium, zero-temperature RFIM spins 
Si = ±1 are placed on a hyper-cubic lattice. At each 
site i a quenched random local field hi is chosen from 
a Gaussian distribution, p{hi) = '^/^^^^P{~ ■ 
field hi acts as a source of disorder for the system, and 
the standard deviation R of the random distribution is 
termed the 'disorder'. The energy of a system with N 
spins is given by ^3 



Jlr \ ^ 



(1) 



where _ff is a tunable external magnetic field. The first 
term of Eq.^couples nearest-neighbor spins ferromagnet- 
ically {{ij) implies summing over nearest- neighbor pairs 
of spins) while the last term provides an infinite range 
anti-ferromagnetic (AF) coupling, where one sums over 
all pairs of spins regardless of their relative distance. 
This AF coupling models the dipolar interactions rele- 
vant in soft ferromagnets, which Zapperi, et. al. 0, 
have shown in three dimensions to have the same effect 
on long length scales as infinite range mean field interac- 
tions. Note that by choosing Jlr = one recovers the 
nearest-neighbor RFIM. In order to model hysteresis we 
study the model at zero temperature, far from thermal 
equilibrium |3 , and for convenience set J = 1 . 

Simulations of the above model were run by starting 
with the external field at H = —00 with all spins down, 
and then adiabatically slowly moving the external field 
through a particular history. As the field H is changed, a 
given spin will fiip (cither upwards or downwards) when 
its effective local field, h"/-^ = H+hi+JJ2{ij) ^j~JLFiM, 
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FIG. 1: Hysteresis loops with concentric inner subloops for 
(a) RFIM with 240^ spins, Jlr = 0.25, and R = 1.8, (b) 
RFIM with 240^ spins, Jlr = 0, and R = 2.7, (c) a Co/Pt 
multilayer thin film, and (d) a CoPtCrB alloy thin film. 



changes sign. Here M = Si is the magnetization of 

a system with iV spins. When a spin flips, it may induce 
its neighboring (or for Jlr ^ even distant) spins to 
flip, creating an avalanche of flipping spins, which is the 
analog of a Barkhausen pulse. The simulations are based 
on the code available on the web 0| which has been 
modified to allow for subloops in the history. The code 
uses the sorted list algorithm, which stores the random 
fields, and is described in depth in Ref. [l^ . 

In the absence of LR fields, as one tunes the system 
disorder (i?), the RFIM exhibits a non-equilibrium sec- 
ond order phase transition at a critical disorder R = Rc 
[^. Below Rc the coupling of nearest-neighbors domi- 
nate the dynamics, and there is a finite jump in magne- 
tization of the saturation loop. However, above Rc the 
random fields dominate, resulting in smooth hysteresis 
curves and mostly small avalanches. Many quantities as- 
sociated with this critical point display scaling behavior 
for R ^ Rc- A detailed discussion is given in Ref. Q. 

The two distinct behaviors of subloops in the RFIM 
can be seen in Figure In Fig.^a) subloops spaced by 
AM = 0.1 are shown for a system with J^n = 0.25 and 
a disorder of i? = 1.8, chosen smaller than Rc = 2.16 as 
many experimental samples such as soft ferromagnets are 
believed to be below the critical disorder. The subloops 
show a linear M{H) behavior as the field is increased, as 
one would expect for a system undergoing domain wall 
propagation while being exposed to a demagnetizing ef- 
fect. Only the first few subloops show some effects of 
the initial condition (all spins down) before the sweeping 
field has created a large enough domain to allow for single 
domain wall propagation. On the other hand one finds 



a very different picture in the absence of LR fields. In 
Fig. n^b), obtained for a system with Jlr = 0, R = 2.7 
and subloop spacing AM = 0.1, as one moves inwards, 
subloops begin to resemble the saturation loop at an 
effectively higher disorder, i.e. R > 2.7. Indeed, pre- 
fiipped spins not participating in a given subloop may 
act as an added, possibly correlated, "effective disorder." 
Tuning the history in this way, we were unable to di- 
rectly observe a transition from loops with a jump in 
magnetization to smooth inner subloops. This is due to 
the inability to break up the system spanning avalanche 
present below the critical disorder even for the largest 
simulated system sizes (480'^ spins) 

Experimental magnetization curves for thin films also 
exhibit these two types of behavior depending on the 
presence of LR fields. Fig. QJc) displays magnetiza- 
tion curves for a Co/Pt-multilayer with the field applied 
along the surface normal. Due to the strong interface 
anisotropy of such multilayer samples, the easy axis of 
magnetization is perpendicular to the film plane even 
though the LR demagnetizing effect is largest in this di- 
rection. Despite the fairly rectangular major loop shape, 
indicating easy axis orientation, the loop exhibits an ex- 
tended linear segment on which all minor loops merge 
due to the presence of LR dipolar effects. On the other 
hand, films with in-plane magnetization behave quite dif- 
ferently, because LR demagnetizing effects are extremely 
small in this geometry. Fig. shows magnetization 
curves for a CoPtCrB-alloy film with the field applied 
in the film plane. This film is polycrystalline with grain 
sizes narrowly distributed around 10 nm diameter and 
exhibits strong exchange coupling between grains. Due 
to the lack of LR coupling, subsequent minor loops ap- 
pear increasingly sheared with decreasing coercive field, 
very similar to the curves shown in Fig. ^b) for the 
Jlr ~ RFIM. All experimental data were measured at 
room temperature using an Alternating Gradient Mag- 
netometer. 

Much information on these two behaviors can be ob- 
tained by examining the Barkhausen noise present in the 
system. With LR fields present, the subloops display a 
power-law avalanche size distribution, all with the same 
exponent and cutoff size, indicating that the system is 
soc B El El- Figure [3 displays the avalanche size 
distribution for subloops of a 21 cm x 1 cm x 30 /im 
ribbon of a Fe2iCo64Bi5 amorphous alloy. In the ex- 
periment, a solenoid provides a triangular driving field 
along the long axis of the sample at 0.03 Oe/s. The first 
cycle drives the sample to saturation, while subsequent 
cycles drive the field to successively smaller amplitudes 
(subloops). The driving rate is slow enough to ensure 
that avalanches do not overlap. The Barkhausen noise 
was measured with a pick-up coil of 150 turns of cop- 
per wire wound around the middle 0.3 cm of the sample. 
By integrating the pick-up voltage one may obtain the 
magnetization. The resulting hysteresis loops are shown 
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FIG. 2: Experimental avalanche size distributions for 
subloops of a soft ferromagnet. The three subloops analyzed 
are shown in the inset, with the largest corresponding to the 
saturation loop. The distributions were extracted from a win- 
dow of width /S.H « 2QA/m which started near H = QA/m. 
A power law of —1.75 is show by the offset, dash-dot line. 
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FIG. 3: Integrated spin-flip correlation functions for 480^ 
systems at _R = 2.198 and averaged over 20 random seeds. 
Curves are given for subloops starting at values of Hmax = 
1.42274, 1.41967, and 1.41317. The inset contains a collapse of 
the three respective distributions, yielding d-\-l3i /ui — 3.0±0.2 
and 1/vf = 1.28 ± 0.40 with H^ax, = 1.427. 



in the inset of Fig. |21 The avalanche distributions for 
the subloops are all identical, affirming the SOC behav- 
ior. The measured power law exponent (f .75) is larger 
than has been predicted and measured for the satura- 
tion loop (1.25-1.5) [ll[l3. This is due to M{H) not 
being strictly linear in the window analyzed. Choosing 
a smaller window close to the linear regime results in a 
power law exponent of 1.3, consistent with previous mea- 
surements. The RFIM with LR fields presents an almost 
identical picture of subloops. For a 240^ system with 
R = 1.8 and J^ji = 0.25, subloops displayed identical 
power laws in their avalanche size distributions with a 
loop integrated exponent of 1.7 ± 0.2, in good agreement 
with the experimental results. 

In the absence of LR fields the case is much different; 
the system history acts as a tuning parameter affecting 
a subloop dependent cutoff. Although a transition like 
the one found in the saturation loop cannot be directly 
observed, its existence may be ascertained by examining 
the scaling behavior on one side of ( "above" ) the critical 
point Ig • Thus we present a scaling collapse of loop 
integrated spin-flip correlation functions and a finite-size 
scaling collapse of the second moments of the integrated 
avalanche size distribution for subloops with J^/j — 0. 

For the spin-flip correlation function, in analogy with 
the saturation loop |^ one may assume a scaling form 

(2) 

where Hmax is the maximum field along the subloop 
and / = Hmaxa — Umax IS the history-induced analog 
of the reduced saturation loop disorder jij]. Here R'j: is 
the critical disorder for the demagnetization curve which 
is numerically within the error bars of Re- Similarly, 
Hmaxa = is the critical field associated with the de- 



magnetization curve [20j. According to scaling theory 
one expects the exponents d + Pi/vi, i^i, and i^f as well 
as the scaling function to be universal. Here a sub- 
script I denotes an exponent associated with the random 
disorder scaling of the critical subloop and subscript / 
one associated with the history-induced "disorder" . For 
the collapses, the system is run at the effective critical 
disorder of the saturation loop, Rc{L), for the linear sys- 
tem size L with Rc{L) — > i?c as L — > oo. This system size 
dependent critical disorder Rc{L) is defined as the disor- 
der at which the maximum number of system spanning 
avalanches in the saturation loop is observed 

The integrated spin-flip correlation function, 
GY*{x,Hmax), was measured for subloops spaced 
by AMmax — 0.025 in their maximum magnetizations 
(see Fig. A similar behavior to the saturation 

loop is found. As the history-induced disorder is 
increased, by moving to inner subloops, the curves 
display a decreasing cutoff. As the system is started 
at the effective critical disorder, Eq. |21 reduces to 
Gf\x,Hmax) ~ x-^'^+f^'/''''^gf\xri). This Scaling 
ansatz results in the collapse shown in the inset of Fig. O 
Only subloops void of spanning avalanches were used in 
the collapses to remove the effects of the finite system 
size, equal to 480"^ spins. 

The integrated avalanche size distribution, 
with a similar scaling ansatz DY^{S, Hmax) ^ 
g-(T+ai3S),jyint^gcrf j^-^^ was prcviously analyzed in 
Ref. 01 . The conclusions therein remain valid, and the 
exponents (see TablePi have been updated to correspond 
with the use of the scaling variable /. 

Table lists the exponents from the above two col- 
lapses along with the saturation loop values. While both 
power law behaviors, given by r -I- af3S and d + P/v, are 
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TABLE I: Universal critical exponents from scaling collapses 
in three dimensions for the history-induced disorder present 
in subloops and the random disorder of the saturation loop. 



Exponent History-Induced Saturation" 

T + apS 2.01 ±0.10 2.03 ±0.03 

1/a 2.3 ±0.5 4.2 ±0.3 

d + f3/iy 3.0 ±0.2 3.07 ± 0.30 

l/f 1.28 ±0.40 0.71 ±0.09 

p 2.6 ± 0.40 2.90 ±0.16 



"References 0,01 

identical within error, the exponents governing the cut- 
off behavior, ct/ and Vf, are found to differ from their 
saturation loop counterparts with only a slight overlap 
in their error estimates. One would expect identical val- 
ues if the pre-flipped spins at the start of a subloop were 
randomly distributed about the lattice, thus preserving 
the uncorrelated, random nature of the system's disor- 
der. However, this is not the avalanches have 
left pockets of unflipped or preflipped spins of all sizes 
up to the system's correlation length. So the difference 
in exponents is not surprising as the history-induced dis- 
order acts more like a correlated disorder as opposed to 
the uncorrelated, random system disorder R. 

The second moments of the integrated avalanche size 
distribution also behave similarly to the saturation loop 
and have a scaling form, in analogy to the saturation 
loop and Eq. 12 {S^)f* ~ L-P'SjifL^/^f) where pi = 
— (t; ± (TiPiSi — 3)/aivi. Table contains the exponents 
from a collapse for systems with sizes from i = 80 to 
L = 480. As each system size was run at the effective 
disorder, Rc{L), a system size dependent Hmax^{L) was 
required in the collapse. The exponent Vf was chosen to 
be consistent across all collapses of the correlation func- 
tion and second moments. Its value differs from the value 
of V for the saturation loop, which reaffirms the differ- 
ences between the history-induced and random disorder. 

Two new exponents, <7f and Vf, were introduced to 
describe the history-induced scaling. However, only one 
is an independent exponent as they obey the exponent 
relation crfi/j: = owi, which may be derived from the 
relations in Ref. [20j. Additionally, the critical subloop 
exponents may be close to those of the saturation loop. 
Indeed, within error bars, numerical results have found 
equal power-law exponents. However, currently no expo- 
nent relations are known that require equality between 
the saturation loop and subloop (subscript I) exponents. 

The RFIM has been used to investigate subloops of 
the main saturation hysteresis loop. Experimental mea- 
surement of subloops on thin film magnets both with and 
without LR forces showed qualitative agreement with the 
RFIM. In the presence of LR forces, subloops may be 
explained by simple domain wall propagation, as con- 



firmed by Barkhausen noise measurements on soft mag- 
nets. However, in the absence of such forces the disor- 
dered critical point was reflected in subloops. Scaling col- 
lapses were performed for integrated avalanche size distri- 
butions, integrated correlation functions, and the second 
moments of the avalanche size distributions. Finally, it 
was shown that only one of the new critical exponents is 
independent of those found in the critical subloop. 
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